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Abstract. The improved data on the cosmic microwave background (CMB) 
anisotropy allow a better determination of the adiabaticity of the primordial 
perturbation. Interestingly, we find that recent CMB data seem to favor a contribution 
of a primordial isocurvature mode where the entropy perturbation is positively 
correlated with the primordial curvature perturbation and has a large spectral index 
(rii so ~ 3). With 4 additional parameters we obtain a better fit to the CMB data by 
Ax 2 = 9.7 compared to an adiabatic model. For this best-fit model the nonadiabatic 
contribution to the CMB temperature variance is 4%. According to a Markov Chain 
Monte Carlo analysis the nonadiabatic contribution is positive at more than 95% C.L. 
The exact C.L. depends somewhat on the choice of priors, and we discuss the effect of 
different priors as well as additional cosmological data. 
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1. Introduction 

A fundamental question in cosmology is the origin of the density perturbation, from 
which the structure of the universe (galaxies, and galaxy clusters) has grown. Today a 
popular scenario is inflation, where the perturbation originates as a quantum fluctuation 
of the inflaton field. 

Clues to the origin of structure can be sought in the nature of this perturbation. 
Simple inflation models produce adiabatic perturbations, but more complicated (multi- 
field) models may also produce nonadiabatic perturbations. 

Current observations are consistent with adiabatic primordial perturbations, but 
sizable deviations from adiabaticity remain allowed. Here "primordial" refers to the 
early radiation-dominated epoch (e.g., at some time soon after big bang nucleosynthesis), 
when all cosmological scales (> 1 Mpc) are well outside the horizon. Adiabatic 
perturbations remain adiabatic while outside the horizon, but give rise to entropy 
perturbations as they enter the horizon. Adiabatic perturbations are completely 
characterized by the associated (comoving gauge) curvature perturbation 71, whereas 
nonadiabatic perturbations have entropy perturbations Sij = + uii) — + Wj) 
between different constituents i, j to the energy density. Here 5 is the dimensionless 
relative density perturbation and w = p/p is the ratio of pressure to energy density. 

A general perturbation can be divided into an adiabatic mode and a number 
of isocurvature modes, which evolve independently. The adiabatic mode has = 
initially (i.e., outside the horizon in the radiation-dominated epoch), whereas 
isocurvature modes have 1Z — and ^ initially. There are four different types 
of isocurvature perturbations [1] - the cold dark matter (CDM), baryon, neutrino 
density, and neutrino velocity isocurvature modes. For simplicity, we consider here only 
the CDM mode, with an initial entropy perturbation 

S = 5 C --^. (1) 

The baryon mode is observationally very similar. 

Even though purely isocurvature perturbations have been ruled out [21 [3] , the data 
allow a subdominant (~10% level [4]) isocurvature contribution. Depending on how the 
primordial 1Z and S perturbations were generated, they may be correlated with each 
other [5]. 

Observationally the CDM isocurvature mode differs from the adiabatic mode in the 
locations of the acoustic peaks in the CMB angular power spectrum Cg. The presence 
of an isocurvature contribution, especially a correlated one, appears in Ce as a change 
in the ratio of peak separation to the first peak position. Since the first peak is well 
fixed by the present data, an isocurvature contribution would appear as a shift in the 
position of the other peaks, reduction in peak separation. CMB measurements 

of increasing accuracy allow thus a tighter constraint on the isocurvature contribution. 
For studies utilizing the Wilkinson Microwave Anisotropy Probe (WMAP) 1-year data 

M , see H m l El HOI ttD EE21 US M EES] • 
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The WMAP 3-year data [TH] is an improvement in this respect [T7J [TBI [TJJ]. While 

the first peak was already measured very accurately in the first year, the determination 
of the second peak shape and location is improved with the 3-year data. Likewise the 
new Boomerang data [20] begins to define the third peak. This motivates a new study 
to update our earlier results 

The focus of this paper is on what the CMB data say about the nature of primordial 
perturbations. Thus we use CMB and large-scale structure (LSS) data only, but address 
other cosmological data in the end. Inclusion of LSS data was needed to break certain 
parameter degeneracies [21] and to constrain extreme values for spectral indices. Current 
CMB data do not cover with good accuracy a sufficient range of scales to constrain 
well the several independent spectral indices of our model. The Planck satellite will 
eventually fix this situation. 



2. Model 



We consider a flat (Qq = 1) ACDM model with primordial curvature and entropy 
perturbations, which may be correlated. For details of our model, see [I]. We give 
below the main points. 

We divide the primordial curvature perturbation into an uncorrelated and a fully 
correlated part. We assume the power spectra of these perturbations and correlations 
can be characterized by power laws, but allow different spectral indices for the entropy 
and curvature perturbation. Thus the spectra can be written as 

V n {k) = Cnn(k) = A 2 r k n ^~ l + A\k n ^-\ 

V s {k) = C ss (k) = B 2 k n -°-\ (2) 
Cns(k) = C sn {k) = A s Bk n ™~\ 

where k = k/ko and ko = O.OlMpc" 1 is the pivot scale at which the amplitudes are 
defined. The spectral index n cor is not an independent parameter and is defined by 
^cor = (^ad2 +^iso)/2. Our approach is phenomenological, looking for signatures of 
deviation from adiabaticity in the observations, and is not tied to a particular model for 
generating the primordial perturbations. The main assumption is that the spectra can 
be approximated by constant spectral indices over the range of distance scales probed by 
the data. For motivation and discussion of our parametrization in the inflation context, 
see [HE21E31E1]. 

The total CMB spectrum Ce can now be divided into four components: the 
uncorrelated adiabatic part, the correlated adiabatic part, the isocurvature part, and 
the correlation between the last two [I]: 

C t = A 2 [(l - a)(l - | 7 |)Cf 1 + (1 - a)| 7 |Cf 2 

+ aCf° + sign( 7 ) v/ a (l-a)|7|Cn ( 3 ) 

/^adl I ^ad2 , ^iso , ^icor 
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where we have defined the total amplitude, the isocurvature fraction, and the correlation 
at the pivot scale 

A 2 = A 2 r + A 2 + B 2 , a = ^ «G[0,1], (4) 
7 = sign(A SJ B)-^-^ 7 6 [-1,1]- (5) 

The Ci denote spectra obtained with unit amplitude (A r = 1, A s = 1, or B = 1). 

Using this parametrization as a basis for a likelihood study has the problem, that 
for small values of a or 7, the spectral indices n; so and become unconstrained. This 
slows down the convergence of the analysis and, if flat priors for the spectral indices are 
used, may bias the results toward zero a or 7. 

Therefore, as suggested in [I], we have switched to using as primary parameters 
the amplitudes at two different scales k\ = 0.002 Mpc" 1 and &2 = 0.05 Mpc -1 , making 
the spectral indices derived parameters. This parametrization and the related issue of 
priors is discussed more in the Appendix. 

In both parameterizations our model has 10 parameters. Four relate to the 
background cosmology: the physical baryon density {uj^ = h 2 Qb), the physical CDM 
density (u c = h 2 Q c ), the sound horizon angle (6) and the optical depth to reionization 
(r). In the old "index" parametrization the six related to perturbations are the logarithm 
of the overall primordial perturbation amplitude (lnA 2 ) and the fractional amplitudes 
a and 7, all taken at k , and the three spectral indices n a( ji, n a d2, and ni S0 . In the 
new "amplitude" parametrization the six primary perturbation parameters are the 
amplitudes In A 2 , In A 2 ,, a±, 02, 71, and 72, taken at ki and ki- 



3. Analysis 

We created 8 Monte Carlo Markov chains (MCMC) using the CosmoMC [25] engine, 
which we have modified to handle correlated adiabatic and isocurvature modes. We 
accumulated a total of 1, 130,306 steps with average multiplicity 10.1, summing up to 
a total of 11, 367, 942 samples. The chains are well mixed, with the Gelman and Rubin 
[26] statistic R — 1 — 0.06. In split (1/2) tests the 95% CL. limits for the primary 
parameters change typically by less than 10% of the standard deviation, and even for 
the worst ones by less than 20%. 

Likelihood of a model was assessed using the WMAP 3-year data and likelihood 
code (version v2p2p2) pH M, EE], Boomerang [20] and ACBAR [29] data, and the 
LSS data from the SDSS data release 4 luminous red galaxy sample (LRG) [301 [21]. To 
avoid the necessity of nonlinear corrections to LSS growth, only the first 14 /c-bands were 
used. We marginalize over the galaxy bias parameter and the ACBAR and Boomerang 
calibration uncertainties. 

Thus in our main analysis we use only four different data sets. The motivation 
for keeping the number of different data sets small, is that with too many data sets 
the interplay of their different systematic effects can become problematic. We use the 
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Table 1. The best-fit models. A: The full model with a correlated isocurvature mode. 
B: The adiabatic model. 
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best available CMB and LSS data only. These are the data where the acoustic peak 
structure can be observed and thus the data which directly probe the adiabaticity of 
the perturbations. So the main object of this study is to see whether the observed peak 
structure is that one which characterizes adiabatic perturbations. 

We assign flat (i.e., uniform over some range; see Appendix) prior probabilities 
to our 10 primary parameters, using the amplitude parametrization. Since some 
parameters are not tightly constrained by the data, the posterior likelihoods depend 
somewhat on the prior probability densities. To show the effect of the priors, and to 
compare to [4], we used the Jacobian of the parameter transformation (see Appendix) 
to convert the 1-d marginalized likelihoods obtained from our chains to correspond to 
flat priors for the index parametrization (dotted lines in Fig. [T|). 

4. Results 

We show 1-d marginalized likelihoods for selected primary and derived parameters 
in Fig. [TJ Instead of the amplitudes at k\ and k 2 we show the derived parameters: 
spectral index rii SO and the amplitude parameters a and 7 at the intermediate scale 
ko = O.OlMpc -1 . The mappings between these are presented in the Appendix. 

The data lead to likelihood peaks at clearly nonzero values for a (the ratio of 
the primordial entropy perturbation power to the total perturbation power at ko) 
and favor a positive 7. (This significantly reduces the pivot-scale dependence of the 
likelihoods discussed in [3].) If the subdominant entropy perturbation is correlated with 
the dominant curvature perturbation it has a much stronger effect on the observables. 
To account for this we also show a COI = sign(7)-\/a:(l — ajpyl, which gives the relative 
"weight" of the correlation spectrum C% OT in Ce, see Equation (jHJ). We find a = 
0.082 ± 0.046 (mean ± stdev), with a < 0.169 at 95% C.L. and a cor = 0.139 ± 0.055, 
with 0.044 < a cor < 0.252 at 95% C.L. The correlation between the isocurvature and 
adiabatic modes is positive, 7 > 0.051 at 95% C.L. 

These values are obtained from the likelihood functions corresponding to the 
amplitude parametrization with flat priors for the primary parameters. As can be 
seen from Fig. [1] the exact confidence levels depend on the assumed priors, but the 




Figure 1. Marginalized likelihood functions for selected primary and derived 
parameters. The solid black curves are our new results using WMAP 3-year data 
and Boomerang data from the 2003 flight (+ ACBAR & SDSS). The dotted black 
curves show the effect of assigning flat priors in the index parametrization instead 
of the amplitude parametrization. The red/gray curves are for an adiabatic model 
using the same data. The dashed blue curves are from our previous study [4] using 
data available in 2004. Note that also in the adiabatic model WMAP3 data favors a 
larger Hq than WMAP1 (not shown) — allowing isocurvature modes favors larger Hq 
regardless of which WMAP data set is used, although the effect is much stronger with 
WMAP3. 



conclusions are not changed at a qualitative level. We discuss the priors more in the 
Appendix. 

Since the definitions of the primary amplitude parameters (e.g. a) depend on the 
choice of pivot scale, we also define 

aT = E(2* + i)c, — ' (6) 

the total nonadiabatic contribution to the CMB temperature variance, 

We find a T = 0.043 ± 0.015, and the whole 95% C.L. range 0.017 < a T < 0.073 is 
positive. Thus the CMB data clearly favor a nonadiabatic contribution. 
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Figure 2. The CMB temperature angular power spectrum for our best-fit model 
(black) compared to the best-fit adiabatic model (red/gray). The dashed blue curve 
shows the nonadiabatic contribution. The inset shows the 2nd and 3rd peaks. 

Our best-fit model and also the best-fit adiabatic model are given in Table I. They 
are compared in Fig. [2j The WMAP 3-year data prefer a slightly narrower 2nd peak 
than the adiabatic model can produce. This holds for the third peak also, but now the 
peak position and width are determined by the Boomerang data. 

This feature in the data can be accounted for by a correlated isocurvature 
component. It can narrow down the 2nd peak without affecting the 1st peak position, 
which is accurately determined by the data. Increasing 9 shifts the whole peak structure 
to the left, making the peaks narrower. Adding a positively correlated isocurvature 
component returns the 1st peak to its place. 

Compared to the adiabatic model, adding 4 parameters improved the fit by 
Ax 2 = 9.7. This comes mainly from the fit to the WMAP and Boomerang temperature 
(TT) Ci\ for the WMAP part from the 2nd peak. The data on the TE and EE Ct are too 
inaccurate for a few per cent isocurvature contribution to play any role. Therefore the 
contribution to x 2 from these cross-correlation data is the same for the adabatic and our 
model. The fit to the SDSS data has improved slightly, whereas that to the ACBAR 
data is indifferent. In Table [2] we give quantitative numbers for the contributions of 
different data sets to the best-fit models. The best-fit model and its x 2 with the data 
are of course not affected by the issue of priors, which affects only the likelihoods. For 



Hints of Isocurvature Perturbations in the Cosmic Microwave Background? 



8 



Table 2. The \ 2 of the fit of the best-fit models to the data, and the contributions 
from the four different data sets. A: The full model with a correlated isocurvature 
mode. B: The adiabatic model. 
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the first time in the history of CMB temperature anisotropy observations, adding the 
CDM (or baryon) isocurvature degree of freedom improves the fit to the data, and the 
likelihoods of isocurvature parameters peak at non-zero values, as seen in Tables [TJ & [2] 
and Fig. [TJ 

The other major effect is that the correlated isocurvature model favors a smaller 
CDM density u c (also Uf, is down, but not as much). This is due to the correlation 
component C^ or , which raises the first and third peaks with respect to the second one. 
See Fig. [2j Lower u> m and u;& compensate this by raising the second peak. 

For fixed ui c and Ub, increasing 9 leads to a larger an d a larger Hubble constant 
Hq. For fixed 9 and cj&, a lower u c requires an even larger Q\ and Hq. Therefore these 
models have a larger H and a smaller matter density parameter Q m = 1 — Q\, than 
the adiabatic model. 

5. Other cosmological data 

The large H = 80 ± 4km/s/Mpc and small Q m = 0.204 ± 0.028 of our nonadiabatic 
model is in some tension with other cosmological data, like the Hubble Space Telescope 
(HST) value for the Hubble constant, Hq = 72 ± 8km/s/Mpc j31j and the estimates of 
n m from Supernova la data (ft m = 0.29+o;al from [32], or Vt m = 0.263 ± 0.042 ± 0.032 
from [33]). To assess the effect of an Q m constraint, we postprocessed our likelihoods by 
weighting (importance sampling [25]) our MCMC chains with corresponding Gaussian 
distributions for Q m . We also ran new MCMC chains with a Gaussian H = 72 ± 
8km/s/Mpc j31j prior. The results are shown in Fig. [3j The effect of the weaker 
Q m = 0.263 ± 0.074 prior is minor and does not eliminate the favoring of nonzero a 
and positive ax- The same is true for the Hq prior, although now the 95% C.L. range 
—0.005 < ax < 0.065 includes a-r = 0. An Q m = 0.30 ± 0.04 prior has a larger effect, 
but the 1-d likelihood of a appears still to peak at a nonzero value. 

As noted in [1], n iso ~ 3 (when n a( j — 1) means that the relative isocurvature 
contribution to Ci is about the same on all scales, but in the LSS power spectrum the 
isocurvature and correlation components overtake the adiabatic ones at small scales. 
Thus there is more (matter) power at small scales, with a larger rms mass fluctuation 
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Figure 3. Marginalized likelihood functions for selected primary and derived 
parameters after including additional priors coming from other cosmological data. The 
solid black curves are the same as in FigfT] The solid blue curves show the effect 
of adding the HST Hq prior. The solid red curves show the effect of the stronger 
fl m = 0.29^0 03 prior and the dashed red curves are for the weaker fl m = 0.263 ± 0.074 
prior. The central values and 1-er ranges of these Gaussian priors are indicated by 
vertical lines and gray areas in the SI a and Hq panels. 



on 8 /i _1 Mpc scale a§. While the adiabatic model favors er 8 « 0.8, our model favors 
as ~ 1. Assuming the power law is a good approximation for the primordial spectra also 
at smaller scales, Lyman-a data can constrain nj SO [31]. Creating new MCMC chains 
with Lyman-a data included like in [34J we checked that the peak in the likelihood 
of rii SO shifted to 2.1, and a to zero. The other likelihoods became more like in the 
adiabatic case. The peak in a T remained clearly positive, but the 68% C.L. region, 
«t = 0.022^q;q25, now included olt = 0. Thus the large small-scale power implied by 
a relatively large constant n iso means that inclusion of small-scale data significantly 
weakens the preference for a nonadiabatic contribution. One possible interpretation is 
that ri; so is not a constant, but becomes smaller for the large k probed by Lyman-a 
data. 
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6. Discussion 

The acoustic peak structure in the latest CMB data favors a few per cent nonadiabatic 
contribution. Our best fit to the CMB + LSS data has a significant contribution from a 
positively correlated isocurvature mode, so that the isocurvature fraction, i.e., the ratio 
of primordial entropy perturbation power to total perturbation power at the scale of 
k = 0.01 Mpc -1 , is 5% and the nonadiabatic contribution to CMB temperature variance 
is 4%. With the addition of four model parameters, the fit to the data improved by 
Ax 2 = 9.7, which can already be considered interesting. 

We note that the nonadiabatic models that give the best fits to the CMB data, 
have some problems with other cosmological data. These problems may be remedied 
if one considers non-flat models. A preliminary look at these models seems to indicate 
that the best-fit non-flat models are closed models with a total energy density a few per 
cent over the critical density, a Hubble constant close to the HST value, Q m and Q\ 
in agreement with the SNIa data, the isocurvature spectral index n; so closer to 2 than 
3, and an even larger nonadiabatic contribution than in the case of flat nonadiabatic 
models. 

While interesting, the preference for nonadiabatic primordial perturbations 
discovered in this study should not be taken as conclusive evidence for a primordial 
isocurvature mode. It is not yet at a very high confidence level, so the possibility 
remains that it is just a statistical fluke. Other possible explanations are an unaccounted 
for systematic effect in the WMAP and Boomerang data, or some other nonstandard 
cosmological feature whose effect on the CMB is in a similar direction. It remains for the 
future more precise CMB measurements to resolve this. If the presence of a primordial 
isocurvature mode were confirmed, it would automatically rule out single-field inflation, 
since generating a primordial entropy perturbation requires more than one degree of 
freedom. 
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Appendix A. Priors and parametrizations 

When some parameters of a model are not sufficiently tightly constrained by the data, 
the posterior likelihood functions become sensitive to the assumed prior probability 
densities for the parameters. Even when one assumes flat, i.e., uniform, priors for the 
primary parameters of the model, the question remains, which parameters are taken to 
be the primary parameters, since the priors for the quantities derived from the primary 
parameters (derived parameters) will not be flat. To avoid problems related to spectral 
indices becoming unconstrained when the corresponding amplitude parameters have 
small values, we introduced in this study a parametrization in terms of amplitudes at 
two different scales, instead of using spectral indices as primary parameters. 

Mapping from the amplitude parametrization to the spectral index parametrization 
is easy to find by employing definitions (j2J), (00), and (jSJ). The spectral indices can be 
written in terms of the parameters of amplitude parametrization as 

, ln[Pgi(fc 2 )/Pgi(fci)] ^ ^^^ 
ln(/c 2 /fci) 

Had2 " 1 " Hhlh) ' ( } 

Hlso 1_ Hhlh) ' (A ' 3) 

where the first (uncorrelated) adiabatic, the second (correlated) adiabatic and the 
isocurvature power spectra at scales ki (i = 1, 2) are given by 

7> wl (A i ) = i4?(l-a i )(l-|7i|) > (A.4) 
Vn2{h) = A 2 (l-ON, (A.5) 
V s (k l ) = A 2 l a i , (A.6) 

respectively. Then the amplitudes A, a, and 7 at ko = 0.01 Mpc -1 are obtained from 
the amplitude-parametrization amplitudes Ai, a±, and 71 defined at fci = 0.002 Mpc^ 1 
by ® 

A 2 = A\ \{\ - a x ){\ - bil)^ 1 " 1 + (1 - ai)|7 1 |jfc n » da - 1 + ajz 71 ^" 1 ] (A.7) 

ai ^"iso-l 

a = = = , , A.8 

(1 - ai)(l - |7i|)A; n adi-i + (1 _ a 1 )\'y 1 \k n >**- 1 + a^^-i 

7= JlH 3 (A.9) 

(1 - |7i|)£; n adi-i + | 7l |^n ad2 -i 

where k = k /ki, and the spectral indices are obtained from Eqs. (1A. ![) — (1A.3[) . 

If we take as a starting point the spectral index parametrization, then the mapping 
from it to the amplitude parametrization is simply 



A 2 = A 2 



(1 -a)(l - \-i\)k naA1 ' 1 + (1 -cOM^ 2 " 1 + ak nis °- 1 , (A.10) 



ak n ' mo 

aii = = = ~ , (A.ll) 

(1 - a)(l - \^f\)k n ^- 1 + (1 - a)\~j\k n ^- 1 + ak 12 ™- 1 
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0/t n ad2 — 1 

7i= — , (A.12) 

n — Jjfc«adl — 1 -)- |<y|fc n ad2 _1 

where fc = ki/ko and i = 1, 2. 

Since we assume that all the component spectra can be described by power laws, 
7i and 72 must have the same sign. Hence, they are not completely independent. To 
obtain independent primary parameters, we draw 71 from the range [—1, 1], but 72 only 
from the range [0, 1], and let 71 determine the sign of the correlation. This has no effect 
on formulas ( 1A.1I) - ( 1A.11I) . but to maintain consistency we need to modify equation 
(1A.12I) . Instead of 7 we should have in the nominator I7I, if % — 2. 

Employing mappings (lA.lj) - (1A.3[) and (1A.7|) - (IA.9I) we obtain the posterior 
likelihoods of n a di, ra a d2> ^iso, A, ck, 7 for a MCMC run in the amplitude parametrization 
(corresponding to flat priors for Ai, a±, 71, A2, 0*2, and 72). Likewise, if a MCMC run was 
made in spectral index parametrization, we could use mappings ( 1A.10I) - ( 1A.12I) to obtain 
the posterior likelihoods of A±, a±, 71, A 2 , a 2 , 72 (and other parameters) corresponding 
to flat prior probabilities for the primary parameters of the index parametrization. 

However, if we want to convert the results obtained in the amplitude 
parametrization to flat priors for the spectral indices, then the mapping flA.lj) - flA.31) 
and ( 1A.7I) - ( 1A.9I) is not enough. We have to correct for the prior too. This can be 
done by weighting the multiplicities in the MCMC chains (i.e. weighting the posterior 
likelihood) by the Jacobian of the transformation ( 1A.1I) - ( 1A.3I) and (1A.7I) - ( 1A.9I) . If 
the MCMC run was made using primary parameters {Oj} (and flat priors for them), 
but we want to show the results with flat priors for {€)«}, the multiplicities must be 
multiplied by 



J 



det 



99, 



(A.13) 

In Figs. IA1I and IA2I we show what the priors of one parametrization become when 
the priors of the other parametrization are taken to be flat. 

In the amplitude parametrization the range of parameters (which is also a crucial 
part of the prior) is throughout our study 

to b G [0.005, 0.1], uj c E [0.01, 0.99], 1000 G [0.3, 10.0], r G [0.01, 0.3] 
ln(10 10 ^) G [1,7], a x e [0,1], 71 e [-1, 1] 
ln(10 10 ^) G [1,7], « 2 G[0,1], 72 € [0,1]. 

In the spectral index parametrization the range of background parameters is the same 
as above, and the perturbation parameters have a uniform prior probability over the 
ranges 

^adi G [-3,4], n ad2 G [-3,4], n iso G [-3, 12], 
ln(10 10 A 2 )G[l,7], «G[0,1], 7 6 [-1,1]- 

These are used only for producing Fig. IA21 since all our analysis presented in this paper 
is based on MCMC runs in the amplitude parametrization. However, we have checked 




Figure Al. The red curves show the priors for the parameters of the index 
parametrization, when flat priors are assumed for the parameters of the amplitude 
parametrization. The black curves are the posterior likelihoods we have obtained (the 
same as the solid black curves in Fig. [1]). 
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Figure A2. The priors for the parameters of the amplitude parametrization, when 
flat priors are assumed for the parameters of the index parametrization. 



against a MCMC run using the index parametrization, that its results agree with the 
dotted black curves in Fig. [I] (which were obtained from the MCMC run using amplitude 
parametrization by correcting for the prior with the Jacobian (1A.13I) ). For Fig. IA2I we 
chose the above ranges of spectral indices to match this check run and our previous 
study [1]. As indices are not symmetric around the scale invariance (n = 1), we can see 
mild asymmetries in Fig. IA2I 

The parameters a cor and «t are derived parameters and therefore their prior 
distributions are not flat in either parametrization. We demonstrate the situation with 
a cor in Fig. lA31 We see that a cor as well as the parameters 77, iso , 7, and a (see Fig. lAll) are 
not well enough determined by the data to make their likelihood functions insensitive 
to the choice of priors. The solid blue curve in Fig. IA3I is the ratio of the posterior 
and prior likelihoods for a cor , and illustrates what the likelihood of a COT could be, if a cov 
had a flat prior probability density The actual likelihood of course depends also on the 
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Figure A3. Prior and posterior likelihoods for the derived parameter a cor . The solid 
red curve is the prior in the amplitude parametrization and the dashed red curve is 
the same in the index parametrization. The solid and dashed black curves are the 
corresponding posterior likelihoods. These are the same as the solid and dotted black 
curves in the a cor panel of Fig. [T] The solid blue curve is the ratio of the posterior and 
prior likelihoods in the amplitude parametrization. 

other priors. 

In Fig. IA4I we show the prior and posterior distributions of ax in the amplitude 
parametrization. We see that the prior is relatively flat in the region where the 
posterior mostly lies. Consequently, unlike for a cor the ratio of posterior and prior 
likelihoods (posterior/prior) for ax is almost indistinguishable from the posterior. Thus 
the situation with a? is better in this respect than the case of a cor , and the conclusions 
about ay are more robust. 
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Figure A4. Prior thin red and posterior thick black likelihoods for the derived 
parameter ax- Both likelihoods are for the amplitude parametrization. Posterior/prior 
would be almost indistinguishable from the posterior. 
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